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INTRODUCTION 


The purpose of this report is to present a simplified computer 
program capable of predicting boiling behavior in aligned liquid face 
seals. The mathematical model used is the same as that used by Hughes 
and Kennedy [1] with the exception that a further simplification is 
made that the vapor region is assumed isothermal. The previous 
results show this to be a reasonable assumption which greatly simpli- 
fies the analysis and calculations. 


ANALYSIS 


This analysis is based on the seal model shown in Figure 1. The 
following assumptions have been made throughout the modeling process: 

1) The liquid flashes instantaneously to a vapor. 

2) The vapor region is isothermal. 

3) The flow is axi symmetric. 

4) Inertia effects in the fluid are neglected. 

5) Heat conduction in the radial direction within the fluid is 
neglected. 

6) The seal plates can be treated as semi-infinite sclids. 

7) The fluid temperature does not vary across the film and is 
equal to tne surface temperature of the seal plates. 


General Equations— Steady State 

The equations of motion which describe the fluid flow between the 
seal plates are 
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for the radial and angular component's. The mass flow rate is 


m = 2n r | pudz = constant 
0 


( 3 ) 


from the continuity equation. The temperature in the film may be 
expressed as: 


-1 / 2 

dT = - r cose) 2 + r 2 sin 2 e] RdRde . 


= ski 


( 4 ) 


which was derived in the previous report, k is the mean thermal 
conductivity of both face plate and nose piece. 

Solving (1) with the boundary conditions that u = 0 at z = 0 and 
z = h yields 


u = kf z(z - h) (5) 

Solving (2) with the boundary conditions that V = 0 at z = 0 and 
V = roi at z = h yields 

V = ru f • (6) 


Substitution of (5) into (3) yields the expression 
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Equation (4) has been integrated for the configuration in 
Figure 2 by Hughes and Osterle [2] and discussed in general terms in 
the previous reports. The results are given below. 


For r-] 1 < r < r 2 ' (liquid region): 

T(r) . ft yi <^) 2 " +4 ] ♦ J o B n [, 


(-~ r ) 2n “ 3 ]) + T 


For r < r-| ' (inside the liquid region) 
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For r > r^' (outside the liquid region): 
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this analysis only the first of these equations is needed. The 
necessaW logical substitutions are presented in Figure 3. 

Assuming the specific volume of a vapor to be much larger than 
that Of ft \iquid and that the vapor acts like an ideal gas, the 


Clapeyron equation can be written 
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The ideal gas law states 


-= RT 


(10) 


Particul a r Equations- - Stead.y State 

The above equations can now be manipulated to define the behavior of 
the liquid. In order to simplify the logic required for inside and outside 
seals two subscripts, i and o, will be used. refers to a dimension 


or property, X, of the seal or fluid entering the seal space, and X 


refers to a dimension or property of the seal or fluid exiting the 
seal space. Thus for an inside seal (P-j > P^) p. = p.j , r. = r-j , P c 
? 2 » r Q = r^, etc. for an outside seal (P-j < P 2 ) P^ = P 2 , r^ = r 2 . 


p o = p 7> r o " V etc - 


Equation (7) can be integrated for the pressure in the liquid 
region (p = constant) 


P - P. = In r/r. 

”4 h 
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or in dimensionless form 


P “P-j ln f/ r -j 
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For the vapor region, assuming isothermal flow at the flash tem- 
perature, the ideal gas equation is used to find p and the integration 
of (7) yields 
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(12a) 


or in dimensionless form 


P 2 - P b 2 _ In r/r b 
P Q 2 - P b 2 ln r o /r b 


(12b) 


P b can now be found by equating the leakage of the vapor in (12a) to 
the leakage of the liquid in (11a) 


( p b - _ < P o 2 - p b 2 > 


u in r,/r. 


W 2 ^g RT b ln r o /r b 


Sol ving 


where 


P b = [A^ + 2AP i + P^] - A 


, . S p g RT b 1n r <Ab 

a " Pin r b /r. 

The Clapeyron equation (9), is used to test if the correct value of 
r b has been chosen. Assuming h^ constant, the integration of (9) 
yields 


Jr^ - exp[- ^ (1 . -!_)] (14) 

p sat K ‘b ‘sat 

where T sab and P gat describe a reference saturation point sufficiently 
close to T b and P bc so that h^ may be considered constant. 
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When P^ is sufficiently close to P^, the leakage and load 
characteristics of the seal may be found. The leakage is found by 
using either equations (11a) or (12a). 

The absolute load, W, supported by the seal is found by 


W = 


b 

r 


P * 2irr dr = 


+ W. 
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(15) 


The total load is simply the liquid load, W £ , added to the load of 
the vapor region, W , P is given by (lib) for the liquid region and 
(12b) for the vapor region. The integration gives 


2 2 (P b ' P i )(r i 2 ' r b 2) 

U. = ,|P b r b 2 - P.r. 2 + 2 in-yrf — 


o o o In r/r. 1/2 
W.. = 2n II [p b 2 + (P 0 2 - P b 2 ) TFT77-] ^ 
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COMPUTER PROGRAM 

logic 

Only the logic of the main program will be presented here. The 
logic of the subprograms is either trivial or stated clearly in the 
previous sections. The program is listed in Appendix A, and it defines 
variable names used. See Figure 4 for the flow chart. 
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Using the Program 

The program is written in Fortran and contains four subprograms: 
(1) A factorial subprogram, (2) A temperature distribution subprogram, 
(3) A Clapeyron pressure subprogram, and (4) A gaseous load sub- 
program. Input is entered via a data deck which is described below. 

The first card of the deck contains fluid information as well as 
plate conductivity. From left to right on this card appears: 

XNAME - the name of the fluid to be sealed (alphabetic) 

MU, (y) - the liquid viscosity (lb-s/ft 2 ) 

MUG, (y g ) - the gas viscosity (lb-s/ft 2 )' 

RHO, (p) - the liquid density (Ibm/ft 3 ) 

XK, (k) - the mean plate thermal conductivity (Btu/hr-ft-°R) 

RVAP, (R) - the ideal gas constant (ft-1 bf/1 bm-°R) 

NHFG - .the number of saturation states to be inputted. 

Figure 5 shows a sample fluid information card. 

The Clapeyron equation, (14), requires that 'saturation information 
be inputted for the particular liquid sealed. Because the Clapeyron 
equation assumes constant h^ gJ the use of only one or two states may 
result in sizeable errors. To correct this, the program reads a satu- 
ration state matrix and uses the saturation data for the point closest 
to T^, the temperature at boiling. There must be NHFG of these points 
supplied (in this case 25 saturation cards were used). The program 
also requires that these cards be placed in order of increasing satu- 
ration pressures. • Each card contains the following three values: 
PSAT(I), (P £at ) - any given saturation pressure (psi a) 

TSAT(I), (T .) - the saturation temperature at the above 

pressure (°F) 

HFG(I), (h f ) - the heat of vaporization at the above 

‘ 9 temperature and pressure (Btu/lbm). 
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Figure 6 shows a sample saturation state deck. 

The final type of data card required is the seal information 
card. From left to right on this card appears: 

TIHF, (T m ) - the bulk temperature (°F) 

PI, (P-j) - the fluid pressure at r-j (psia) 

P2, (P 2 ) - the fluid pressure at r 2 (psia) 

R1 , (r-j) - the smaller seal radius (in) 

R2, (r 2 ) - the larger seal radius (in) 

H, (h) - the thickness of the lubricating film (y in) 

OM, (u) - the angular velocity of the seal (rpm) 

There is no limit to the number of seal information cards which can 
be entered at one time. The only requirement is that the final card 
must have TINF (TJ = 10,000. Figure 7 shows a sample seal information 
deck. 

The program uses the data and outputs pressure and temperature 
distributions, leakage rates, and absolute load. The program is 
capable of handling liquid seals and gas seals in addition to mixed- 
phase seals. All output is given in both English and SI units. The 
computer program and sample output is found in Appendices A and B. 


RESULTS AND CONCLUSIONS 

The computer program (Appendix A) was used to analyze two seal 
configurations. The first configuration approximates that used by 
Orcutt [3] (r-j = 2.025 inches (0.0514 m) ; r 2 = 2.225 inches (0.0565 m)). 
The second configuration approximates a commercial seal manufactured 
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by the Crane Packing Company of Morton Grove, Illinois (r^ = 1.693 
inches (0.0430 m); r 2 = 1.849 inches (0.0470 m)). The thermal con- 
ductivity, k, is chosen to be the average of the thermal conductivities 
of the two seal plates. For example, in the case of the Orcutt seal 
k was taken to be 7.5 8tu/hr-ft 2 -F (13.0 W/m-K), (kq Uar j. 2 * 4.0 
Btu/hr-ft Z -F (6.9 W/m-K); fe carbon . graph1te i H.O 8tu/hr-ft 2 -F 
(19.0 W/m-K}}. To find the effect of k on loading, other values were 
also used. The results of the computations are summarized in 
Figures 8-19. 

The following conclusions have been made: 

1) When phase change is included the pressure distribution is 
radically different from the simple linear pressure distri- 
bution commonly assumed in industry. ■ The forces actually 
pushing the plates apart is greater due to phase change 
(Figures 8 and 9). 

2) Leakage decreases when flashing occurs. The leakage, from an 
inside seal is approximately the same as the leakage from an 
outside seal (Figure 10). 

3) The model breaks down when boiling occurs close to the 
entrance. This apparently happens because the model assumes 
instantaneous flashing, but when saturation conditions are 
reached very shortly after entrance into the seal, an 
extended region of mixed vapor-liquid flow probably exists. 

This is in agreement with Orcutt' s observations (Figure 11). 

4} For a given seal configuration there are two film thicknesses, 
h, which support the same load, W. When the spacing under- 
goes a small excursion about the larger equilibrium value (due 
to a vibration or some other external disturbance), the seal 
will return to its original position. However, an excursion 
about the lower space is not stable. An increase in spacing 
grows until the larger equilibrium spacing is reached, but a 
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decrease in spacing causes a catastrophic collapse of the 
seal faces. Metal -to-metal contact or rapid explosive 
boiling might occur as the spacing decreases. The con- 
sequent explosion then forces the plates apart. These 
phenomena have been observed under certain conditions 
(Figures 12, 13 and 14) . 

5) For a given seal configuration (and given spacing) there are 
two angular velocities, w, which produce the same load, W. 

The seal is stable at the upper to but unstable at the lower to 
(Figures 15 and 16). 

6) For a given seal configuration there are two bulk temperatures, 
T^, which produce the same load, W. The seal is stable at the 
upper T m , but unstable at the lower T m (Figures 17 and 18). 

7) Varying k changes the boiling radius and shifts the load 
curves either left or right. Plates with a high k will con- 
duct more heat from the fluid and vaporization will occur 
closer to the exit. A low k has the opposite effect. 

8) Equation (8) shows the seal load to be a’ function of the 

2 2 
parameter h/cu . There are two ratios of h/w which support 

a given seal load. One is stable and the other is unstable. 

(Figure 19). 
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NOMENCLATURE 



m 

e 

p 

P 



V 


w 

W 

Z 


- distance separating the seal plates 

- heat of vaporization 

- thermal conductivity of the seal plates . 

- mass flow rate 

- circumferential coordinate 

- density 

- pressure 

- pressure found by the Clapeyron equation, (14) 

- rad.ial coordinate 

- ideal gas constant of the vapor 

- temperature 

- radial flow velocity component 

- absolute viscosity of the liquid 

- absolute viscosity of the vapor 

- circumferential flow velocity component 

- angular velocity 

- the absolute load supported by the seal 

- axial coordinate 


Subscripts 

1 - a dimension or property of the seal or fluid at the inner 

radius 

2 - a dimension or property of the seal or fluid at the outer 

radius 

b - a dimension or property of the seal or fluid at the boiling 
radius (liquid-vapor ir.-rface) 
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1 - a dimension or property of the seal or fluid at the point 

where fluid enters the seal space 

z - pertaining to the liquid region 

o - a dimension or property of the seal or fluid at tfVd point 
where fluid leaves the seal space 

sat - defining a saturation state 

v - pertaining to the vapor region 

00 - pertaining to the bulk properties 
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Hughes and Qsterle Model for Temperature Distribution on 
the Surface of a Semi-Infinite Solid. 
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Figure 3. Logical Substitutions for the Temperature Distribution 
Equation, . (8a). 



Figure 4. Program Flow Chart. 
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Figure 5. Fluid Information Card for Water. 



Figure 6. Saturation State Deck for Water. 
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Figure 7. Sample Seal Information Deck. 
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Figure 8. Effect of Phase Change on the Pressure Distribution. 


i 



PRESSURE, 


3C9 

(422) 


(103) 1.700 
(4.3 1 8) 

Exit 


1.740 

(4.420) 


290 

1.780 1.820 (4 16) 

(4.521) (4.623) 

Inlet 


RADIAL POSITION, r. inches (m x \ 0 ? ) 


Figure 9. Temperature and Pressure Distribution in the Film Space. 

P ) = 15.0 psia (103.4 kPa), P 2 = 90.0 psia (620.4 kPa), 
ri = 1.693 inches (0.04300 m). r 2 = 1.849 inches (0.04696 m) , 
h * 50.0 x 10" 6 inches (1.27 x 10' G m) , u> = 7200 rpm 
(754 rad/s), k = 26 Btu/hr-ft-°F (45 W/m-K). 
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Effect of Bui k Temperature on Leakage. = 45.0 psia 
(310.2 kPa), P 2 = 15.0 psia (103.4 kPa), rj = 2.025 inches 
(0.05144 m), r 2 = 2.225 inches (0.05652 m) , h = 50.0 x 
10 -6 inches (1.27 x 10” 6 m) , u> = 3600 rpm (377 rad/s), 
k B 26 Btu/hr -ft-°F (45 W/m-K). 
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Figure 11. Fluid Temperature and Pressure at Various Seal Spacings. 

Pj = 15.0 psia (103.4 kPa), P 2 = 90.0 psia (620.4 kPa), 
ri = 1.693 inches (0.04300 m) , r 2 = 1.849 inches 
(0.04696 m), T„ = 270°F (405 K) , u = 7200 rpm (754 rad/s), 
k = 26 Btu/hr-ft-°F (45 W/m-K). 
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SEAL LOAD 




ro 

O'* 



FILM THICKNESS, h, inches x I0 6 (m x I0 6 ) 

Figure 13. Seal Load at Various Film Thicknesses for an Outside Seal. P, = 15.0 psia (103.4 kPa) 

x 2 = oS;^ P /! a ^ 6 ? 0 * 4 kPa )‘ r i = 1 - 693 inches (0.04300 m), r 2 = 1.849 inches (0.04696 m), 
T. = 270 ° F ( 405 K), a) = 7200 rpm (754 rad/s), k = 60.2 Btu/hr-ft-°F (104 W/m-K). 






Figure 14. Seal Load at Various Film Thicknesses for an Inside Seal. 

P! = 45.0 psia (310.2 kPa), P, = 15.0 psia (103.4 kPa), 
n = 2.025 inches (0.05144 m) , r 2 = 2.225 inches (0.05652 m) , 
T = 230°F (383 K) , u = 5000 rpm (524 rad/s), k = 60.2 
Btu/hr-ft-°F (104 W/m-K). 
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Figure 15. Seal Load at Various Speeds for an Inside Seal. P, = 45.0 psia (310 2 kPa) P 
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BULK TEMPERATURE. T.. °F (K) 

Figure 17. Effect of Bulk Temperature on Seal Load for an Outside 
Seal. Pj = 15.0 psia (103.4 kPa), P 2 = 45.0 psia 
(310.2 kPa), rj = 2.025 inches (0.05144 m) , r 2 = 2.225 
inches (0 05652 m), h = 50.0 x 10* 6 inches (1.27 x 10" 6 m) , 
u * 3600 rpm (377 rad/s), k = 26 Btu/hr-ft-°F (45 W/m-K). 
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Figure 18. Effect of Bulk Temperature on Seal Load for an Inside Seal 
Pi = 45.0 psia (310.2 kPa), P 2 = 15.0 psia (103.4 kPa), r, 
2.025 inches (0.05144 m) , r 2 = 2.225 inches (0.05652 m), 
h = 50.0 x 10 -6 inches ( 1 . 27 x 10 -6 m), u = 3600 rpm 
(377 rad/s), k = 26 Btu/hr-ft-OF (45 W/m-K). 
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COMPUTER PROGRAM 
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THE MAIN PRoGrAM STArjs HERE 

THIS PROGRAM SOLVES THE ALlQ N ED RADIAL Face seal PRoBLE m ASSU- 
MING the -VaPOR REgIONI IF ANY>|S ISOTHERMAL. 

THROUGHOUT This Program Variables STARTING WITH The letter R,P, or 

I RE_FEH_JT-0_R AD I i-» _PRLSSU_ReS,_OR_.TE«PERaTURES_aT_V^R J OUS-i’OINTS 

along the SEAL PLATES. 

THE POINTS are DEFINED BT THE FOLLOWING FIVE SUbScRjPTS; 

1 - A VALUE AT ThE INNER RADIUS 

2 - A VALUE AT ThE OUTER RADIUS . - 

B - A VALUE AT THE BOILING RADIUS 

i - a valu e at the point where fluid enters the seal plates _ [ 

o - a value at the point where flu?d leaves The seal plates. 

DIMENSION p I 2 1 ) , To I 2 | ) »RM ( 2 I » . Tqm ( 2 1 ) . Pm ( 2 1 ) 

COMMON/aREA I /R I »RO#RB *R2 ,R 1 .NhAX , AN ( 1 oo ) »BN( IOO ) ,XMU »OM , XK ,H |T 1 NF/ 
IaREA2/PsATI50)»TSaT(S0).HfG(S0»»RVaP#NHFG/aREa3/R(2i),DRG,c1»C2 
Fahrenheit to kelvin 

FT0CIa»»( A-32. >/l ,e>273* 

PSI TO KPASCAlS 

PTOPA ( A ) *A*&.893 ' 

Inches To METERS 
XIT0CMI A I?A*2.SM/100. 

LIQUID pressure DISTRIBUTION 

__ PLCA»«PiMPB-Pi_)*.(ALOG.(A/RI I/aLOG(RB/RI ) > 

BOILING PRESSURE aT INTERFACE 

PIF(A)=SQRT(IEZ*TB*AL0G(R0/Al/A , .OG(A/Rl)>**2.»2.*EZ»Te«Pl*AL0G(RO/ 

IA)/ALOG(A/rI)^Po*POI-IEZ»TB*ALOG(Ro/A)/ALOG(A/RI)) 

GaSeOUS PRESSURE DISTRIBUTION 

pg(ai=sqRT{Pb*pb-»(po*po-pb*pb)*(alog(a/rb)/alog(Ro/rb))) 

MASS_-FLQJ*_ I N _JHE L I QUj D E G 1 0 N 

XMLlA.B)»« I“P1 E*Rh0*h**3. J/(6.«XMU*aL 0GCB/RI » ) »•< A-PI )/l00 0 000.**3 

. mass flow In the vapor region. . . .. . 

XMV(AtB) a -(IP0 # P0-A*Al*PlE*H*»3.)/(i2.*XMUG*RVAP*TB*AL0GlR0 / B)J/10 

100000. *.3 ... _ . . . ..... 

LOAD OP the liquid region 

* L I A . B I » A B S I P I E • I a • 6 • B^_P I _«R_I_! R I ♦ ( ( a -P 1 J • .< RJ ^Rj “ B * Bj J i < 2 . l 0 G < B / R I 

1 ) ) ) ) 

I NPUT_: ... ... 

XnamE - the name of the fluid to be sealed 


• r '''i- .h i wu n; i jq v. / • • 

of The vapor (Lbf-sec/ft**2 i 


i i / t i d a i i j » » n » vj ' l l 

make .input oata dimensionally correct ..and _supplt_const ants 
XK=»CXk«778 . / 3600 • 

PIE*3,IhIS9 ... 

NMAX-IOO* 

_£P.M AX* ,2S 
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7* - 

C 

»* 

C 

H 

c__ 


C 


C - 

2 * 


J* 

• . . . 



»• 


6 • 


;• 


8* 



C 

0 * 


l» 



l» 


3* 


1* 


5* 

12 

»• 


7* _ 


4* 

16 

»• 

20 


• * HERE : _ 

p I e is The mathematical constant 

NMax_JS_THE NuMttE« or TERMS JO_£E-Ja*EN_I 

EPM*X |S THE MAXIMUM eR r0r BETWEEN THE CL 

PReSSUReS(LBE/FT*»2) - 

XMUM*XMU*97 »883 

XMUGM»XHUG*M7.0e3__ . . --- 

RHOM»Rho*1A»018 

RVaP-M = R V AP*S .38 

CX(CM»CXk*I .7303 

RRlTE(A»900)XNAMEtXMU,XMUM»XMUG.XMU6M,RH0 
WRI TEC A »905 ) Ry AP .KVAP m * CXK t CXKM 
SOLVE FOR temp. CONSTANTS ( An ( I ) and bniii 
00 20 J = 1 , nMAX 

X J*Jr_J 

IFIJ .GE. I 2 ) GO TO 12 

F?1FACTI2« # XJ)*«2./FACT(XJ)**‘I*I/2. ••(*♦. *XJ 

GO TO 1 A 

E“<( 8 ./(pIE*l 2 »*XJ*l.)J)*t<XJ*l.)«ExPll.)/I 
l j* | •)/(?«*xJ+2* ) • 1 1 • i l • ♦ 1 • / ( I 2 • • I 

. 21j/(12.«(Xj^1 •III.. 

ANC J)»F/ I 2*x J*H. ) 

BN( J)rF/(2*XJ-3. ) 

Input variable input dataivjd): 
ctinf - the bulk temperature ( f ) 

CPI - THE F|_U | D PRESSURE AT The INNeR RAO 

' CP2._=„THE_FlU ID.PKESSURE. aT_JTHe OUTER. Rad 

cri - the inner radiusiiN) 

- . -CR2 - THE OUTeR RADIUSUN) 

CH - The FLUID ThICKNeSSCmJcROjNCHES) 
com - the Seal angular veloc i t y ( rpm j 

REA0I5, n20)CTINF. CPI |CP2,CR1»CR2, CH.COM 

IF l A BS ( C-T- UJ F 1 0 0 0 0 • J- • L T 0 0 1 J_S T Op 

EZ=RHO*xMUG*RVAP/XMU 

CONVERT VID To PROPER UNITS, __ 

T 1 NF«CT ! NF*HAO • 

Pl i -CPL«1*H..__ 

P2*CP2* 1 . 

R4 “XR4-/-I-2# 

R2-CR2/1 2. 

M*CH/4-2-» 

OM*C0m*2»*pIE /& 0» 

5EJ —UP -SEAL- AS -INFLOW OR -OUTFLOW 

IFIP2 .gT. PI >G0 TO 3 S 


N THE T 
aPetron 


EMP£R 
A Nq 


ature 

I TERA1 


. Rhom 


2 • • X J ♦ 1 
2 • • X J ♦ 1 


. J ) • 
• ) > ) 


1 2 . • I I 

< • 2 • ) / 


2 • • X 
< I .♦ 


I US ( PS l 
1US(PSI 


A ) 

A) 


PI*P1 

R0-R2 

P0»P2 

GO _T 0 _R 0 

R 1 -R2 

PI «P2 

RU«R I 

P0-P1 

TEST FoR a complet 

R B “ R 0 

T0-TI90) 

1FJ PBC-I-j OJ -PnJ 6 0 • AO J.4S 


complete liquid SEAL ( I TTpE- I I 
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I 


12* 


0 
J* 

4 » 

if. 

»• 

|7* 

J* 

»• 

3* 

If. 

2 * 

)• 

1* 

IS* 

I** 

7* 

3* 

?• 

0* 

I* 

2* 

3* 

1* 

5* 

3* 

7* 

8* 

9» 

0* 

I* 

2 * 

3* 

1* 

15*_ 

6 • 
70 
30 

I**.. 

0* 

I* 


. c 


ISO 


160 


• ■*HPB,R8|*WGlPB,RBl 
. . . CONVERT TO PROPER OUTPUT UNJTS 
CT 1 NFM»r TOC I CT INF ) 

_jCPJM«PT 0 PA(CP1J 

CP2M-PTOP A ( CP 2 ) 

crim*xitocm«cri ) 

CR2M«X I T OC M ( C R 2 ) 

CHM-XIToCM(CH) ... 

C0HH«c0m*2.*PIE/60» 

— RB =RB *4 2 

tb-tb-i&o. 

-pb*pb/ir**» - - * — 

00 ISO 1*1.21 

. _ R C I )»R ( I ) ‘12. - - 

p I I ) - P < I ) / 1 M H • 

„ _ R M (1 ) * X 1 T OC H < R < 1 ) ) 

PH C J ) »PToPA 1 P ( I ) » 

XMM»XM*.MSM 

fth=ft*M .H*<a 

•rite the Pertinent information 

»R1TEU»910)CT1NF*CTINFM.cP1 ,CP!M,CP2,cP2M»cR 1 ,cR 1 M ,CR2 

..•R|TEU»9»S)Ch*ChN.C0m,C0mM 

irllTYPC »eQ» l ) WR I Te * 6 »?20 > 

1F1ITYPF .E0« 2)«r1TE‘6,92S> 

IFUTYPf »E0* 3)i(lrRlTE<6f9307 

IFUTYPF • E 0 • 2 ) Go TO 170 

• R I TE ( 6 »93S ) 

IRU£16»?JLU 

00 160 1=1,10 

•RITE(6i9H6)R(I),RM»1),TD(I),Tdm(I>.,P(I),PM(I) 

IFUTyPe • E 0 • I 1 WR 1 TE 1 6 19S6 ) R « 1 1 ) , rh ( 1 1 ) , TD (1 1 1 , TqM ( 1 1) ,P ( j 

IFUTYPe • EO • 3)«rRlTE<6»9S0)R< 1 1 ) »RM« 1 1 1 ,TD«1 1 » ,T D M1 l 1 ) ,P( i 

_U_) 


, C R2H 


1 I . P« C 1 
1 ) ,PM( 1 


1 70 


ISO 

190 


IF ( ITyPe .eQ» 1)60 TO 190 

__ »RIJ£(6 ,9SS) 

• R 1 T E ( 6 »9M0) 

DO 180 1=12.21 


2* 
3 • 

c 

S. 

c 
_ c 

&• 

7* 

800 
8 1 0 

!• 

fj» 

820 

900 

0* 
1 • 
2 • 
3* 

90S 

1* 


S* 


6* 

910 

7* 


8* 


9* 



GO 


•RlTE<6.99S)RCl),RMll),P<j7,PM(I) 

•RItE I 6 i960) X M , X m Nj ft_, jf, M 

TO 3c 


FORMAT STATEMENTS 


FORMAT(a1S,2E10.3,3F10.3.iS) 
E_Q RMAT L3£J-0iiD 


FORMAT ( 7F1Q.3) 

FORMATJy ! JHE_ FlU.I 0 _T0 BE S E A L E D Is i_l l A 1 0/ / / , •_ MU . ThE.LJqUID VI 

ICOSITY m • , E 1 0 • 3 . * LB _ S/Ft** 2 *».E10.3. f PA^S •//• MuG. THE GA 

2 _V I SC OS I T Y ■ » .E I 0. 3 .» LB-5/F T •* 2 . a !_, E 1 0 • 3.* *_ PA-S * // • rHO » THE 

3 LIQUID DENSITY = •,F10.3,' LfiM/FT • • 3 ■ »,F10.3,* K(/M«»3«/l 

FORM AlX! R.V.AJL* _LHel_J D L A L_g A S_C Q N S T A N T_ * Fl . 2 , A_F_T rL 8 F y LB m. -R - = 

1 F 7 • 2 . ' .i/KG-K'//' THE THERMAL CONDUCTIVITY OF THE SEAL PLATE * ' , F 

27.2.! .BtU/hR-FT-F = !_J .FT .2 . • ft / M - k _•_) 

FORMATI'I T1NF * »,F6.1 t * Dt G F = •,F6.|.' 01 
I • PSI A.* • ,F7. I , • KPA _P2_- . i .F6. 1 . • RSI* 


2 * *»F6.3.' 

3 !J 


G K * / * P| ■ » , F 6 • 1 , 
• » F7 • 1 , K p a ' ✓ ' R 1 


I N 


' , F 7 , S 


M 


R2 


iF6.3,« IN = • , F 7 . 5 • # 


OC 


V. 


ti w i 






= R I 

»RB 

= 0 • 

»0 » 

BS ( a'QR 
_20 L« l 
XL-L-1 

>•__ AA« AN t LJLv? 1 I O R I / A 1 * *!2 • «XLM .Jj 

I I 20 SUHAsslIHA* AA 

• 1 F C A B S ( a - Q r 2 ) »LT . _«DOOO 1 ) G0_J_0 J.IJO 

• 1 1 30 00 1150 L = 1 .NMAx 

■± X L- L- ) 

• BB=BNIL)««1.-(A/QR2)**(2.*XL-3.I) 

• 1150 SUMB s SL»MB'*BB 

• I 1 *»U T=(110CC > 0On* # XMU # OM*OH # A»*3l/(2 # *XK*H)l* t SUMA^SUMB)>TlNF 

• END 





gffK o»UOW!';« Y u J 0( l “ l 
ORIGIN M 1 Ajh 
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C ” CLAPETrOn PRESSURE SUBPROGRAM 

FUNCTION PbC< A I . 

COMMON /aREa2/PSaT(S0),TSAT(S0J,HFG(S0)»RVAP*NHFG 

1FIA . L T « (TSATtl)+H& 0*l) PBC*PSA T(M *EXP< ( - H F G ( 1 _L* 27 8 » / R V A P 1* < 1 ./A . 

1»1*/CTSaT( I 1 ♦H 60 • ) ) 1 • I *• M • 

1 F ( A < GT « IT SAT { NhFG I ♦.M60. )_) PB Cji PS A T ( NHFG ) *E X P ( ( - HF 6 < NHFG I *7 7 8 . /R V 

1aP)*<1./A-|./(TSAT<NhFG1-M<sO»)))MmH. 

I F ( < A .1 T. ( TSA T ( 1 ) ♦R60 » ) ) .OR« (A . GT A .< T 5 A T ( N hF G ) a H fc 0 . ) ) ) R E T UR N 

L-2 

1200 I F ( c A .pE« (TSATI l IH6Q.)) *A nD« ( a « Ll»_! I S XL ItULL* M 6 0 - 111 GO TO 
1210 

L »L*1 _ . 

GO TO 1200 

12 1 0 L1«L 

|FU A-HaO*-TSATIL1 > .GT. 1TSAHL^1)-(A“*<60.)11L1=L*1 

PBt c PSAT (LI ) * E X P ( (-HFGIL1 ) • 77 & » /R V AP ) • < 1 « / A- 1 » / ( T S A T I L 1 ) ♦ M & 0 » > ) > • 1 

1 *♦ M . 

RETURN _ 

END 




ORIGINAL PAGE IS <- - 
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G A S E 0 U S L 0 A D_5 U B p R 0 G R A M 

FUNCTION WG*A,B) 

C O Mb OtU.A RE A 3 /-R -L2-1-J -» D R G_. C1_^C 2 

DIFFERENTIAL load of the GASEOUS region 

F»G < A ) !LS O R XL* JC l.tLC 2* A L 0 &J A/RB >JL?A 

R B * B 

Z»F*GIR(l' U.±£Jt6i« L24.1JL 

DO 1 300 L= I 2 ,20 ,2 

130U ( R_l L ! J 

DO 1310 L=l3,20,2 

1310 Z«Zt.2.*_FwG ( R CL-U 

«VG*ABSI?.<3»l‘HB9*DRG*Z/3.) 

RETURN 

END 


Appendix B 
SAMPLE OUTPUT 


A . M 

E FLUID TO BE SEALED IS: 

Water ... _ 

• the liquid viscosity « 

.590-05 LB-5/FT • • 2 

6 » t h e gas viscosity = 

•252-06 LB-S/FT*«2 = 

3» THE LIQUID DENSITY » 

A0»000 LbM/Ft**3 * 

\P. the IDLaL G A S CONSTANT 

« 85.70 FT-LbF/LbM 

: THERMAL conductivity of 

thf seal plate => 7 
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•283-03 PA-S 


• I 2 l -O m PA-S 
961 .080 X6/M • *3 
-R = H 6 \ .07 J/KG-K 

i 

♦ | 

.50 BJU/HR-FT-F =» 12.98ft/M-X 

I 

I 



> 


TINp » 205.0 DEG F • 369.7 DEG K 

PI ■ **5.0 PS I 6 ■ 310.2 K P A P2 > 15»0 PS 1 A * 103. H *PA 

Rl * 2*025 IN ■ .05113 M R2 ■ 2.225 In » *05651 M 

H ■ 50* 0 MlCROlNCHES ■ 1.27 MIcRONS 

OH ■ 1000.0 RPM . 1U1.72 RAD/SEC - . . 


THIS SEAL ACTS AS A LIQUID SEaL 


the liquid distribution i s ; 


R, IN(M) 


T ,F (K ) 


P.PSIA(ICPA) 


2.0251 .05113) 209.35(371.53) 15.001 310«18) 

_ 2.015C. .0519H) 209.55(371 .61 ) _.\ . . 1 1 . 8 7 ( .. 2 8 8 . 6 1 ) 

2 . 06 5 ( .05215) 209.60(37 1.7 1 ) 38.77 ( 267.21 ) 

2.085 ( .05296 ) 209.76 (371.76) 3S.70( 216*08 ) 

2.1051 .05317) 209.80(371.78) 32.661 225.12) 

2.12SC . 05397 ) 209*8 2(371. 79) 29. .65 ( 201 • 36j _ 

2 • 1 *4 5 ( • 05 M 1 8 ) 209.80(37 1.78 ) 26.66( 183*79 ) 

2 . 1 65 ( .05199) 20.9«_7 5 (371 .25) 2 3 .7.1 ( 1 63.11 > 


2.1851 .05550) 209.65(371.69) 20.78< 113*22) 

2..205J * 056 0JJ 209.50(37.1 . 6 J ) 17-.88.1. _1 23 . 22 i _ _ 

2.225 ( «05<>5 1 ) 209 . 26 ( 3 7 1 . 1 8 ) 15.001 103.39 ) 


the Leakage rate = . 177-01 lbm/sec = .802-05 kg/sec 

id E_ A BlS O i. U_I e_L 0 AO .S U PP OR.IE.D_ R y —T H I S_S E A L-. » 7 9 L B F. » 


- 351 • N 



TINf ■ 205.0 OEli f ■ 369. 

PI ■ **5.0 PS I * - JlO.2 KPA 

R l ■ 2 • 0?5 IN ■ .051*0 M 

deg k 

P2 • 15.0 PSIA 

R2 • 2.225 In - 

■ 1 03 . *4 

•05651 

KPA 

M 

H ■ 

OM « 

SQ.O MlCROlNCHES ■ 1.27 MI c RONS 

_ 5000*0 RPM s 523.60 RAD/SEC 

1 


this seal acts as a mixed-phase seal 

... . - - - 



ThE LIQUID DISTRIBUTION'S 

• 




• 





R • I N ( H ) 

T,F|K) 

P ,PS I A 1 KPA ) 

• 

2.025( .05 1 *4 3 ) 

2 • 03 7 C . 05 1 7 h ) 

2 7 3 . *4 2 1 *4 0 7 . 1 2 ) 
_2?6 • 1 6 ( *408 .6*4 1 

*45.00 1 
*4 *4 . 8 8 1 

310*16) 

309.32) 

A 

2 • OH 9 ( .05205) 

2*061 ( .05236) 

277.991*409.66 ) 
_279.18<*410.32) 

*4*4 .75 1 
. _*4*4 .63 1 
*4*4.51 1 
*4 *4 .J) 8 ( 

308 .*47 ) 
307.62 ) . 
306.77 ) 
305.93 ) 

• 

2.0731 .05266) 

2 • OB 5 l .05297) 

279.831 <1 10.68) 
280*01 1*410. 78) 


2.098( .05328) 

2.1101 .05359) 

2 79.72 4*4 1 0.62 ) 
2l 8 . 9 3 1 <1 1 0 . 1 8 ) 

*4*4 .26 1 
*4*4 . 1 *4 ( 

305. 10) 
30*4 • 27 ) 


2.1221 .05389) 

2 • 1 3 H 1 « 05 *4 2Q ) 

277.581*409. *43) 
275.51 1 * 4 0 8 • 2 8 ) 

*4*4.02 ( 
*43.90 1 

303 • *4*4 ) 
302 * 62 ) 

KB = 

2.1*46 IN « .05*451 M 

IB = 272. *4 *4 r = 

*406 . 5b If 


PB «= 

*43.78 PSIA * 301.81 KPA 



the 

Vapor distribution is: 




• 


R, IN(M) 

p ,PSI A 1 KPA ) 




2. 15*4 1 .05*47 i ) 

2.J621 .05*49 U 

*4 1.77 1 287.95 ) 

^39.6 7 1 273* *4S) 

1 



2.1701 .0551 1 ) 

2.1781 . OS 5 3 1 ) 

37. *46 1 258.1 9 ) 

35.111 7*42.03 1 




2.185 1 .0555 l ) 

2.1931 .0557 1 ) 

3 2.6 1 ( 22*4 * 78 ) 

^29.9 1 ( ^206 .16) 




2.201 ( .0559j ) 
2. 2091 .0561 i ) 

26.951 185*76) 

2 3 • 6 3 ( 1 62*91 ) 




2.21 7 1 .0563) ) 

2.22 5 1 .0565 1 ) 

19.791 136.38) 

15.001 103.39 ) 



the 

leakage rate » . 116-05 lbm/sec = . 528-06 

KG/SEC 


the 

absolute load supported by 

This seal » los. 

lbf « 

*465 . N 


BEPR0DUC1BIU1Y UK iM 
OBIGiNAT PACE IS P<X* 
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TlNf = 2bO • 0 DEG ► « H|n ffl |)[Gn" 

P 1 » H 5.0 PS|A s JlQ.2 KP« p 2 = 

K1 = 2*025 IN = .CS1M3 M R2 = 

H » Sq.O MICRO I NCHES * 1.27 MORONS 

°.*1_?__5.Q00.Q_RPM_«_ _5.2.3.6Q_RAD/5£C 

I 


Thlb seal acts as A gas seal 


15.0 PS 1 A 
2*225 In * 


L0 3>H .KPA- 

• 056 S 1 M 


THE V A PO K 0 I S t R I BUI I UN is: 


R , IN ( M ) 


p l PS I A ( K P A ) 


2 • QH 5 ( 
2.0651 

2 . 005 ( 
_2 «_L0 5J _ 
2.1251 
2 • 1 H 5 ( 
2.1651 
-2. 1 85 I_ 
2.205 { 


.05 1 9s ) 
•052MS ) 

.05296 ) 
.OS 4 JLZJ- 
. 0 5 J 9 7 ) 
.OSHMQ ). 
. 05^99 ) 
,GS5So )- 
.05601 ) 
.05 6544- 


h2 .86 ( 2 9 5 » H 5 ) 

-HO . 63 { 2 bo • 101. 

38 • 30 ( 2 6 M • 0 1 ) 

-J5.8JU 2-*l7 .£54_ 


33.22 ( 
30 .-Hi t 
27 .3H ( 


229 -0 1 1 
209 .62 )- 
1 88 . H5 ) 


2.1851 ,QS55o ) 23 « 9 1 ( i 6 q . 8 1 ) 

2 . 205 { .05601) 1 9 • 9 H ( 137.HH) 

2,225 1 — .056.543 j _5 . G 0 ( 1-03 *39 1 

the leakage rate = _ .H 70-r b ‘ lbm/sec = . 213-06 kg/sec 

THE A b SQ L UTE LOAD SUPPORTED B y TuIS SEAL = 86. L B F = 
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